Orientational order and alignment of elongated particles induced by shear 
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Shear induced alignment of elongated particles is studied experimentally and numerically. We 
show that shear alignment of ensembles of macroscopic particles is comparable even on a quantitative 
level to simple molecular systems, despite the completely different types of particle interactions. We 
demonstrate that for dry elongated grains the preferred orientation forms a small angle with the 
streamlines, independent of shear rate across three decades. For a given particle shape, this angle 
decreases with increasing aspect ratio of the particles. The shear-induced alignment results in a 
considerable reduction of the effective friction of the granular material. 



Flow of large ensembles of elongated objects - often 
observed in nature or industry - usually induces pro- 
nounced alignment of the building blocks. This phe- 
nomenon is found at all length scales, in log jams on 
rivers, in slowly deformed multicomponent clastic rocks, 
in seeds, nanorods, viruses, and even at molecular scales 
in nematic liquid crystals [H-Q3- n one hand, such 
alignment processes are poorly characterized for macro- 
scopic objects, even though granular flows have been ex- 
tensively studied in the last two decades [l4]-[29| . On the 
other hand, shear alignment and collective reorientation 
dynamics is well documented and exploited at molecular 
scale (for nematic liquid crystals) and can be described 
quantitatively by continuum theory [30|,|3l|. Studies have 
also shown that the stationary alignment observed for 
simple one component systems may be replaced by a so 
called tumbling behavior for more complex materials, e.g. 
nematic liquid crystals with pre-transitional smectic fluc- 
tuations l3ll - l33| , polymeric liquid crystals (33437| as 
well as dilute suspensions where individual particles ro- 
tate in shear flows (Jeffery orbits, [38]). For the limiting 
case of very dense suspensions direct interparticle forces 
dominate, leading to stationary flow alignment, analo- 
gously to what is expected for simple one component 
systems. 

In this study we focus on simple one component sys- 
tems and aim to bridge the gap between systems of very 
different sizes and interaction types. Our goal is to deter- 
mine the shear induced order and average orientational 
angle for macroscopic objects, to track how the phe- 
nomenon depends on the particle aspect ratio, and find 
to what extent the applied shear rate is relevant. Start- 
ing from an initially random orientation, we also aim to 
quantify the change in the shear resistance of a sample 
due to the evolving alignment. 

We sheared the granular material in a cylindrical split 
bottom geometry (see Fig. [1]) where stationary flow is 



obtained by a continuous rotation of a thin plate be- 
low the material. The strain is localized [26|, [27], [39| in 
the so called shear zone (indicated with dark gray). In 
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FIG. 1: (a) Schematic illustration of the split bottom shear 
cell [26[] . The granular material consisting of elongated grains 
with mean diameter d and length L is sheared due to the 
rotation of the thin circular plate below the material with 
diameter D p « 80 d and angular velocity Q. The container 
diameter is D c ~ 140 d, while typical filling height is H « 
lid ~ 0.14 Dp (unless stated otherwise). Shear deformation 
is localized at the region marked with dark gray, (b) Particles 
are oriented on average at a small angle with respect to the 
streamlines, (c) Photographs of the particles studied, with 
L/d = 2.0, 3.4, 4.5 (rice), as well as 3.5 (glass cylinders) and 
3.3, 5.0 (pegs), (d) Measured borders of the shear zone for 
pegs with L/d = 5.0 (border defined at approximately 10% 
of the maximum shear rate). 

this geometry the time and ensemble averaged shear rate 
changes with the radius in a relatively wide region allow- 
ing us to map the shear rate dependence of the align- 
ment in a single experiment [28|. Prolate particles with 
different length to diameter ratios L/d were investigated 
(Fig. Ht). Two methods were used to visualize and eval- 
uate shear-induced alignment: (i) complete 3D imaging 
of the sample with X-ray tomography and (ii) optical 
imaging of the particles at the upper surface of the sys- 
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tern with a digital camera. For the characterization of 
changes in effective friction due to the alignment, the 
torque on the bottom plate was measured at stationary 
rotation. The experiments were complemented with nu- 
merical simulations using the Discrete Element Method 
(DEM). In the simulations elongated particles were pre- 
pared by numerous overlapping spheres and the align- 
ment of the particles was measured in plane Couette flow. 

The shear-induced orientational order is monitored by 
diagonalizing the symmetric traceless order tensor T: 
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where is a unit vector along the long axis of particle 
n, and the sum is over all N detected particles. The 
largest eigenvalue of T is taken as the primary order 
parameter S. A second, biaxial order parameter D is 
defined as the difference of the two other eigenvalues of 
T. The shear-induced alignment is characterized by the 
average alignment angle # av measured with respect to the 
shear flow direction (see Fig. []})), i.e. the deviation of the 
largest main axis of T from the streamlines. 

First, we present the shear rate dependence of the av- 
erage alignment angle # av and the order parameter S ob- 
tained by digital optical imaging at the surface layer. 
In this experiment, the shearzone at the top was divided 
into 15 bands (widths 0.19 cm « d) and the average shear 
rate was determined for each band. The orientation and 
in-plane length of colored tracer particles (6% of the sam- 
ple) were detected with a particle tracking code, yielding 
about 60,000 particle positions in each band. The center 
# av of the orientational distribution was determined by 
fitting a Gaussian, and the order parameter S was cal- 
culated. These two quantities are shown for the roughly 
ellipsoidal rice grains in Figs. [2^i,b as a function of the 
shear rate 7. Here the shear rate spans a wide range of 
inertial numbers, about 0.0002 to 0.4, using 3 different 
rotation rates of the bottom plate (ft = 0.26, 0.52 and 
4.2 rad/s) and considering the spatial variation of 7 in 
the shear zone. The ranges of 7 spanned for each ft are 
indicated with horizontal bars in Fig. [2^. (The inertial 
number (2l| is defined as / = jdy/ p/P, where p is the 
density of particles and P is the pressure.) 

As it is seen, av is well defined, and it is essentially 
independent of the shear rate for the the whole range of 
7 investigated here. This result is not unexpected, it is 
actually similar to the case of flow aligning nematic liq- 
uid crystals consisting of elongated molecules, where the 
alignment angle is a material parameter, i.e. shear rate 
independent as well. But even if the aspect ratios of the 
building blocks (molecules and granulates) are similar, 
the interparticle forces have very different character in 
these two systems. While the molecules experience at- 
tractive forces (dipole-dipole, van der Waals) the grains 
investigated here interact only via hard core repulsion 
and friction. The fact that despite of these differences 
the alignment angle is independent of the shear rate un- 



derlines the robustness and the geometrical origin of the 
shear alignment. The slightly decreasing value of the 
order parameter S with increasing shear rate (Fig. (2Jd) 
shows that the inertial effects make the material less or- 
dered while preserving the same alignment direction. 

The numerical simulations were performed using the 
discrete element method (DEM), where long particles 
were approximated by chains of overlapping spheres. 
Special care was taken to reduce the possible interlocking 
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FIG. 2: (color online), (a) The average orientational angle 
# av of the particles with respect to the streamlines and (b) 
the order parameter S as a function of the shear rate 7. Both 
panels correspond to optical measurements for rice grains with 
length to width ratios L/d = 2.0, 3.4, and 4.5. The value of 
# av is roughly constant, and S varies only very weakly across 
three decades of the shear rate. 

of neighboring chains: both the overlap (30-70%) and the 
number of participating spheres were varied randomly. 
For implementing a soft particle model we used the gen- 
eral purpose molecular dynamics code LAMMPS |40| , 
extended with particle cementing interactions within the 
chains [29| without allowing the contacts to break. An 
assembly of 600-1200 elongated particles was confined un- 
der constant pressure between two parallel, rough walls 
of finite mass; the system was periodic in the other two 
directions. The shear flow was realized by imposing a 
fixed, opposite velocity of the walls. In the steady state 
the shear velocity profiles were reasonably linear to ne- 
glect shear rate variations. The inertial number was in 
the range 0.004 to 0.09. For calculating the average align- 
ment angle and order parameter only the central 60% was 
considered, in order to exclude influences of the walls. 

The complete 3D arrangement of all particles in the 
shear zone was also determined experimentally for two 
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samples (pegs with L/d = 3.3 and 5.0) using X-ray to- 
mography [46|. The purpose of these experiments is to 
show that the order and alignment data extracted from 
surface particles are representative for the global order in 
the shear zone, with only slight variations. The particles 
were detected by applying a watershed algorithm to the 
recorded images. The components of T were determined 
for a surface region of ~ 1 cm (two particle diameters) 
depth and for a 2 cm thick layer below (bulk), corre- 
sponding to about 10,000 and 20,000 particle positions, 
respectively. 

Figure [3] combines the bulk and surface results from 
X-ray tomography, the surface data from optical experi- 
ments, and the numerical results. Comparing the average 
values of the shear alignment angle for the six samples 
investigated, we find a systematic decrease of # av with 
increasing length to width ratio (Fig. [3^l) for a given ma- 
terial (rice or cylinders). At comparable aspect ratios, 
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FIG. 3: (color online). The average orientational angle # av as 
a function of the (a) length to diameter ratio L/d of the par- 
ticles and (b) order parameter of the system obtained by ex- 
periments (#,4, A) and numerical simulations (x). One datum 
point (a) was taken with a particularly deep layer, H/d = 34. 
The alignment angles for typical flow aligning nematic liquid 
crystals MBBA and 5CB (*,<>,+) @4lil4l| and for a nematic 
liquid crystal 8CB (□) which shows tumbling when approach- 
ing the smectic A phase transition are shown for comparison. 



the alignment angle appears to be smaller for cylinders 
(as if they had a larger effective aspect ratio). Fig. [3b 
shows that the alignment angle # av displays a decreasing 
trend with increasing order parameter both for experi- 
ment and simulation. The measurement with H/d = 34 
shows that the value of av is independent of H. The 
orientational order slightly increases with decreasing H, 



this may be an interesting subject for further investiga- 
tions. As a possible explanation for the trends observed 
in Figs.[3K,b, we suggest the following picture: in a dense 
granular system, the elongated particles are entangled 
with neighbors that hinder their rotation. In a shear flow 
an orientation distribution develops due to the competi- 
tion between the agitations by collisions and this hinder- 
ing effect. The orientation angles of the particles have 
a distribution centered around # av , where the width of 
the distribution roughly equals # av . For more elongated 
and therefore more entangled particles the angular fluc- 
tuations will be smaller, so the width of the distribution 
is smaller, making # av smaller. 

The X-ray tomograph measurements confirm that the 
shear alignment is very similar at the surface and in the 
bulk, with slightly larger order parameter S and slightly 
smaller alignment angle av in the bulk (see Figs. [3ji and 
b) compared to the surface. The slightly larger order ob- 
served in the bulk is connected to a smaller standard de- 
viation from the mean orientation within the shear plane 
(while the out-of-plane standard deviation is similar in 
the bulk and at the surface). This is in good agreement 
with molecular dynamics simulations of a hard-spheroid 
fluid under shear flow [42] , where a similar biaxiality was 
detected. For our case the measured biaxiality is some- 
what stronger in the bulk (D « 0.09) than at the surface 
(D«0.03). 

For comparison, we have also included the order pa- 
rameter dependence of the alignment angle for typical 
flow aligning nematic liquid crystals (MBBA and 5CB) 
for which the aspect ratio of the molecular shape is es- 
timated to be L/d = 2.5 11]. For molecular nematics 
the order parameter is independent of the shear rate (in 
this case determined by the temperature). We emphasize 
that the decreasing tendency of the av (S) curves is very 
similar for granulates and for these nematics. We also 
note that different flow behavior is observed for some 
nematics that exhibit strong smectic fluctuations when 
approaching the smectic A phase transition [l(J l3ll - [33| . 
For example, nematic 8CB shear aligns only in a nar- 
row temperature range near the clearing point. With 
decreasing temperature, the 6 av curve approaches zero 
when one of the viscosity coefficients changes sign, 5.6 
K above the transition to smectic A (where S « 0.4). 
Pretransitional smectic clusters drive the system into a 
tumbling behaviour, where the director continuously ro- 
tates in the shear flow. Another example for tumbling is 
the case of polymeric liquid crystals, where the presence 
of poly meric chains strongly modifies the flow behavior 
[3J-[37|. For molecular systems the flow behavior of the 
aligned case as well as the dynamical modes of the tum- 
bling state can be accurately characterized by modern ex- 
perimental tools, such as X-ray scattering |3J-[36|, NMR 
|34i |37| or conoscopy [32| , and non-equilibrium molecu- 
lar dynamics simulations [33, |42j . Various theories have 
also predicted the instabilities and different dynamical 
modes, like tumbling, kayaking, wagging, log rolling or 
flow aligning (431445). These theories are more general in 



the sense that they map the phase diagram not only as a 
function of aspect ratio but also as a function of density 
including the case of dilute systems. T hey predict flow 
alignment for our parameter range [43|, |44]]. Altogether, 
the fact that in our dense granular system, with no de- 
tectable smectic fluctuations and simple interactions be- 
tween the particles stationary alignment is observed, is 
in agreement with the above observations and theories. 

The initial evolution of order has been measured start- 
ing from a randomly oriented configuration. Data were 
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FIG. 4: Evolution of the orientational distribution of the 
particles starting from an initially random orientation. 

obtained from 1300 experiments. The orientational dis- 
tribution is presented in Fig. 2] as a function of the strain 
taken in the middle part of the shear zone for a sample 
with L/d = 3.5. As it is seen, order develops quickly and 
the alignment angle # av decreases towards the stationary 
value. The strong decrease of # av at the beginning of the 
process is a consequence of the fact that the vast majority 
of particles rotate in the same direction due to the shear 
flow. The transient in Fig. [4] evidences that an equilib- 
rium alignment is established after the sheared particles 
have passed approximately three to four neighbors. 

The increasing orientational order and the develop- 
ment of the shear alignment angle leads to a systematic 
change of the effective friction of the sample at the begin- 
ning of the process. This can be quantified by measur- 
ing the torque M needed to rotate the circular bottom 
plate with a constant speed. As a reference we have also 
recorded the torque for nearly spherical mustard seeds. 
This helps us to correct for the transient effects which 
are not directly related with orientational changes. The 
signals were corrected by subtracting the torque needed 
for rotating the empty plate and rescaled by the sta- 
tionary value Moo. The rescaled torque M = M/M^ 
is shown in Fig. [5^, while Fig. [5]o presents the curves 
for the elongated grains after subtraction of the data 
for mustard seeds on a log-lin scale. The evolution of 
the internal friction of the material can be character- 
ized by an exponential decay of the excess torque as 
M — Mmustard ~ (M max - l) e -^— . The solid lines in 
Fig. [SJd are exponential fits, they define a characteristic 
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FIG. 5: (color online), (a)-(b) Evolution of the applied 
torque needed to maintain stationary rotation of the bottom 
plate starting from an initially random orientation. Panel (a) 
also includes a reference dataset obtained for spherical parti- 
cles which is then subtracted from the other datasets result- 
ing in the curves presented in panel (b). (c) Characteristic 
strain 7trans corresponding to the transient and (d) ratio of 
the torques M ma x/Moo = M max needed to rotate an unori- 
ented sample and the oriented sample (stationary state) for 
rice (•), and cylinders (a). 



deformation 7t ra ns describing the transient. The depen- 
dence of 7trans on the aspect ratio is shown in Fig. [5fc. 
For rice grains, the transients increase with grain length. 
This evidences a slower dynamics of the more elongated 
particles. A linear fit yields 7trans = 0A(L/d - 1) for 
rice. The ratio of the torques corresponding to the totally 
random and the oriented systems is estimated around 
M max & 1.5 as seen in Fig. [5]1. 

In conclusion, we used surface optical imaging, bulk 
X-ray tomography, as well as numerical simulations to 
study the shear flow of elongated dry granular particles. 
We found that the particles quickly orient themselves 
around a preferred direction, which is at an angle from 
the streamlines. This angle is almost independent of the 
shear rate (or inertial number) across three decades, and 
for a given class of shapes (rice or cylinders) it decreases 
with increasing particle aspect ratio. This shear align- 
ment is a robust phenomenon, and its geometric origin 
is underlined by the similarity with nematic liquid crys- 
tals, despite the completely different interparticle forces. 
The shear alignment decreases the effective friction by 
about one third when compared to the initial randomly 
oriented state. 
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